WMAP Constraints on Quintessence 
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We use recent results from the Wilkinson Microwave Anisotropy Probe (WMAP) for the locations 
of peaks and troughs of the Cosmic Microwave Background (CMB) power spectrum, together with 
constraints from large-scale structure, to study a quintessence model in which the pure exponential 
potential is modified by a polynomial factor. We find that the model is compatible with all the 
recent data for a wide range of cosmological and potential parameters. Moreover, quintessence is 
favoured compared to ACDM for n, » 1 and relatively high values of the average fraction of dark 
energy before last scattering ("early quintessence"); for n s < 1, quintessence and ACDM give similar 
results, except for high values of early quintessence, in which case ACDM is favoured. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Recent cosmological observations suggest that there is 
a dark energy component to the energy density of the uni- 
verse, which should be added to the matter component 
so as to reach the critical density. Theorists have consid- 
ered various possibilities for the nature of this dark en- 
ergy, notably a cosmological constant and quintessence, a 
dynamical scalar field leading to a time- varying equation 
of state parameter, = Pa/PiP- These models most 
often involve a single field H @, H S S E BJlEJl3 
or, in some cases, two coupled scalar fields [ill fl3j . 
Other possibilities for the origin of dark energy include 
the generalized Chaplygin gas proposal [PJ] and cardas- 
sian models [P5| . 

In order to unravel the nature of dark energy, it is cru- 
cial to use observations so as to be able to discriminate 
among different models. In particular, the existence of a 
dark energy component affects the structure of the CMB 
power spectrum, which is particularly sensitive to the 
amount of dark energy at different epochs in cosmology. 
For instance, the locations of peaks and troughs depend 
crucially on the amount of dark energy today and at last 
scattering as well as the dark energy time-averaged equa- 
tion of state, which are model-dependent quantities |16| . 
Hence, one can use the high precision measurements re- 
cently obtained by the BOOMERanG [HI, MAXIMA-1 
[l8| . Archeops [l!j and, in particular, WMAP [gjj obser- 
vations to constrain dark energy models. 

This study has already been performed for the case of a 
cosmological constant [2l| , the generalized Chaplygin gas 
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[22j and for some of the quintessence models that have 
been used in the literature, e.g. for the pure exponential 
potential [TEEl, a "leaping kinetic term" model OHl], 
Ratra-Peebles potential [II 111, a class of SUGRA 
potentials [2^ and cosine- type quintessence j2|J. 

The goal of this paper is to study the effect of a dark 
energy component defined by the quintessence potential 
proposed in Ref. [27} 
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on the location of the first three peaks and the first trough 
of the CMB power spectrum (for a first study of CMB 
anisotropies for this model see Ref. [2^ ) . We also analyze 
the consequences of cluster abundance constraints on as, 
the rms mass fluctuation on scales of 8 h~ 1 Mpc. 

The M-theory motivated potential of Eq. leads to 
a model with some interesting features, namely there are 
two types of attractor solutions giving rise to an accel- 
erating universe today, corresponding to permanent or 
transient acceleration |29| . Transient vacuum accelera- 
tion is a particularly appealing scenario that would also 
solve the apparent incompatibility between an eternally 
accelerating universe and string theory, at least in its 
present formulation, given that string asymptotic states 
are inconsistent with spacetimes that exhibit event hori- 
zons [30|. 

For both types of solutions, there is scaling of the 
densities early in the expansion history, with w$ — 
1/3, fi<£ » 4/A 2 in the radiation dominated era and 
Wcj, = 0, rifj, w 3/A 2 in the matter dominated epoch, fol- 
lowed by vacuum domination and accelerated expansion. 
Transient vacuum dominated solutions arise for AX 2 > 1 , 
in which case the potential has no local minimum or, for 
AX 2 < 1, if the field <fi arrives at the local minimum with 
enough kinetic energy to roll over the potential barrier 
and continue rolling down the potential. Permanent vac- 
uum domination occurs for AX 2 < 1, if </> gets trapped in 
the local potential minimum. 

Our study of the dependence of the first three peaks 
and first trough locations on the potential parameters 
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in the (Cl m ,h, n s ) cosmological parameters space, in 
view of WMAP's results, reveals that, for n s w 1, the 
quintessence model of Eq. JU is favoured as compared 
to the ACDM model, provided A, which determines the 
average fraction of dark energy before last scattering, sat- 
isfies A < 15. For n s < 1 and A < 15, we find the opposite 
i.e. ACDM is favoured as compared to quintessence. As 
A increases, A > 18, the model's results become com- 
parable to ACDM's, independently of n s , as should be 
expected since the average fraction of dark energy be- 
fore last scattering becomes negligible. Moreover, the 
model presents a non-negligible fraction of dark energy 
at last scattering and during structure formation, typi- 
cal of early quintessence models, leading to suppressed 
clustering power on small length scales as suggested by 
WMAP / CMB /large scale structure combined data jU . 



II. LOCATIONS OF PEAKS AND TROUGHS 

We consider a spatially-flat Friedmann-Robertson- 
Walker (FRW) Universe containing a perfect fluid with 
barotropic equation of state p w = wp w , where it; is a 
constant (w — 1/3 for radiation and w = for dust), 
together with a scalar (quintessence) field with potential 
given by Eq. Q . The evolution equations for a spatially- 
flat FRW model with Hubble parameter H = a/a are 



H = 


1 (A 
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subject to the Friedmann constraint 
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+ p m + -<j> 2 + V 
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We work in units where Mp = (&-kG)~ 1 / 2 = % = c = 1. 

The CMB peaks arise from acoustic oscillations of the 
primeval plasma just before the Universe becomes trans- 
parent. The angular momentum scale of the oscillations 
is set by the acoustic scale £a which for a flat Universe 
is given by [j^ 



Ia = tt 



TQ — T\s 
C s Tls 



(4) 



where and r = J a" 1 dt is the conformal time, To and t\$ 
being, respectively, the conformal time today and at last 
scattering; c s is the average sound speed before decou- 
pling 



C, = T, 1 



c s dr 



q , 9 Pb (t) 



(5) 



with pb/Pf the ratio of baryon to photon energy density. 

In an idealised model of the primeval plasma, there is 
a simple relation between the location of the m-th peak 
and the acoustic scale, namely £ m w ml a- However, 
the location of the peaks is slightly shifted by driving 
effects and this can be compensated by parameterising 
the location of the m-th peak, £ m , as in |34Ll35| 



i Pm = lA{m- tpm) =t A (m-(p- Sip m ) 



(6) 



where (p = tpi is the overall peak shift and S(p m = ip m — <p 
is the relative shift of the m-th peak relative to the first. 
Eq. JBJ is correct also for the location of troughs if we set, 
for instance, m = 3/2 for the first trough and m = 5/2 
for the second trough. Although is not in general possi- 
ble to derive analytically a relationship between the cos- 
mological parameters and the peak shifts, one can use 
fitting formulae that describe their dependence on these 
parameters. Doran and Lilley |35| give accurate analytic 
approximations for the first three peaks and first trough, 
which can be found in the Appendix for convenience. No- 
tice that, as the authors of Ref. [35j point out, although 
these formulae were obtained using a standard exponen- 
tial potential, one expects the results to be fairly inde- 
pendent of the form of the potential (since the shifts are 
almost independent of post recombination physics) un- 
less it is qualitatively very different from the exponential 
potential before last scattering. 

The locations of the first two acoustic peaks and the 
first trough, from the WMAP measurements of the CMB 
temperature angular power spectrum are [2(j 



£ P1 = 220.1 ±0.8 
£ P2 = 546 ± 10 , 
£ dl = 411.7 ±3.5 



(7) 



notice that all uncertainties are la and include calibra- 
tion and beam errors. The location of the third peak is 
given by BOOMERanG measurements [3(| 



£ P3 — 825^^3 



(8) 



We have studied the location of the first three peaks 
and first trough in the (f2 m , h, n s ) cosmological param- 
eter space for the potential Eq. , for different values of 
A and A, slicing through the (O m , h) an (f2 m ,n s ) planes. 
For each value of A and A, <fio is chosen such that fltot — 1 ■ 
Gi 

[2G| . we iissuiix' i KTcn i [■<•]■ 

Our analysis in the (f2 m , h) plane shows that, for n s — 
1, the strongest constraints arise from the positions of the 
first and second peaks and, for clarity, we plot only the 
£ Pl , £ P2 contours corresponding to the WMAP bounds on 
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FIG. 1: Contour plots of the first and second Doppler peaks location in the (£l m ,h) plane for the potential Eq. Q, with 
A=0.0025, for different values of A and spectral index n s — 1. Full and dashed contours correspond to WMAP's bounds on, 
respectively, £ P1 and £ P2 , see Eq. 0. The dotted contours correspond to erg constraints. The intersection of all curves leads to 
the shaded allowed regions. We also show the corresponding contours for the ACDM model with n s = 1; the full and dashed 
boxes indicate, bounds from WMAP only and a combination of WMAP and other experiments on h and Q m h 2 , Eqs. @ and 
IIUI I respectively. 



these quantities, Eq. Q), for different values of A and A = 
0.0025 (see Fig. [I); also shown are the £ Pl ,£ P2 contours 
corresponding to the ACDM model for the same values 
of the cosmological parameters, where the dashed box 
represents WMAP's bounds on h and Q m h 2 , namely [2(j 



n m h 2 = 0.14 ±0.02, h = 0.72 ±0.05 



(9) 



III. CONSTRAINTS FROM a s 

We have also studied constraints resulting from erg , the 
rms density fluctuations averaged over 8/i _1 Mpc spheres. 

Ref. [33 gives an estimate of the CMB-normalized erg- 
value for a very general class of quintessence models from 
Q$(a), w ( j ) (a) and the o^-value of the ACDM model with 
the same amount of dark energy today = fi9(A): 



The full box corresponds to the bounds obtained on these 
quantities from a combination of WMAP data with other 
CMB experiments (ACBAR and CBJQ, 2dFGRS mea- 
surements and Lyman a forest data [20j 



fl m h 2 = 0.135 



+0.008 
-0.009' 



h - 71 + 04 



(10) 



— ~ (a eq ) f (1 - n A 



N-(l + ffi- 1 )/5 



(11) 



where a eq is the scale factor at matter-radiation equality, 
£l s / is an average value for the fraction of dark energy 
during structure formation i.e. during the matter domi- 
nated era, before starts growing rapidly at scale factor 

O-tr 



For n s — 0.95, the strongest constraints come from 
WMAP's bounds on £ Pl ,£ P2 and £ dl , see Fig.0 

Finally, the results of our analysis in the (fl m ,n s ) 
plane, with h = 0.71, are shown in Fig. |3| where we 
plot the same contours as for Fig. 



An a tT 

fl s J = [lna tr — lnaeg]^ 1 / ^^(a) eHna . (12 

J In a eQ 



The effective equation of state w is an average value for 
W0 during the time in which Q.^ is growing rapidly: 
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FIG. 2: As for Fig. 1 but with n s — 0.95 and the contours corresponding to the observational bounds on the first trough, 
Eq. 0, are also shown (dot-dashed contours). 



f ii^(a)w 1 (a)dlna 

-1 J ln a tv 

— HI • 

I £i,p(a)dlna 

•J In a+ r 



(13) 



In order to compute as for the ACDM model, we use 
the definition 



(14) 



with r = 8 h 1 Mpc and 



A 2 (fc) 



H 



3+n 



T 2 (fc), 



(15) 



where T(k) is the transfer function describing the pro- 
cessing of the initial fluctuations, for which we use the 
fitting function [3^ 



T(q) 



ln(l+2.34q) 
2.34q 

[l + 3.89q+(16.1q) 2 + 
(5.46g) 3 + (6.71 9 ) 4 ] _1/4 



(16) 



with \3£ 



/!f Mpc 



r = Q m h exp 



-O ft 1 



2h 



- - r, 



(17) 



and Sh is the density perturbation at horizon crossing. 
A fit to the four-year COBE data gives [4(j 



10 5 ^ = 1.94a 



~0.785-.051nO„ 



exp [an + bh 2 



(18) 



where n = n s — 1, a — —0.95 and b = —0.169 (assuming 
there are no gravitational waves). Taking into account 
rcionization effects, Eq. 118|l should be corrected; we use 
the fitting formula of Ref. [llj 



5h(t) 
Sh(t = 0) 



= 1 + 0.76r - 1.96t^ + 1.46t j 



(19) 



where r is the optical depth. This formula is reliable up 
toTfs 0.5; we use r = 0.11, which is within the range of 
WMAP's bound, r = 0.166±g;^f HJ. 

We compare our results for as with large scale struc- 
ture data. The recent study of the mass function 
of 300 clusters at redshifts 0.1 < z < 0.2 u sing 
early SDSS data yields a 8 Q° m 6 = 0.33 ± 0.03 
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FIG. 3: Contour plots of the locations of the first and second peaks and first trough in the (O m , n) plane for the potential 
Eq. with A=0.0025, for different values of A and h — 0.71. Full, dashed and dot-dashed contours correspond to WMAP's 
bounds on £ P1 , £ P2 and £d 1 , Eq. £J, respectively. The dotted contours correspond to as constraints. We also show the equivalent 
plot for the ACDM model with h = 0.71. 



Other cluster analysis (using eg. BCS, REFLEX and 
ROSAT data) yield different values: er 8 = 0.43fi m 38 

[13,(78 = (0.508 ± 0.019)fi- (0 - 253±0 024) 0, a 8 = 
(0.7 ± 0.06)(fi m /0.35)-°- 44 (r/0.2) 08 [H and a 8 = 
0.38fi-°- 48 +°- 27n ™ HI. 

The results of our study of a 8 constraints for the model 
of Eq. are shown in Figs. 1 — 3 (dotted curves), where 
we have plotted the largest erg contours that are compat- 
ible both with the large scale structure fits mentioned 
above and the CMB computation, Eqs. l(TT|l - (|T§1) . 



IV. DISCUSSION AND CONCLUSIONS 

Fig. n shows that, as A increases, the allowed region 
becomes smaller and more similar to the ACDM model 
results; in particular, for A = 15, the allowed region is 
only marginally larger than the one for ACDM model, 
which is quite small. 

Notice that a lower bound on A already exists from 
standard Big Bang Nucleosynthesis (BBN) and the ob- 
served abundances of primordial nuclides, which implies 
f^(MeV) < 0.045, or, considering a possible underesti- 
mation of systematic errors, the more conservative result 
f^(MeV) < 0.09 0; these bounds require, respectively, 
A > 9 and A > 6.5, for the model we are considering. 



For n s = 0.95, Fig. shows that the allowed region 
becomes larger until A ~ 15 but does not change sig- 
nificantly for 10 < A < 18. Again, we find that, as A 
increases, A > 18, the model's results become more sim- 
ilar to the ones for the ACDM model. 

Regarding the results of our analysis in the (fl m ,n s ) 
plane, with h — 0.71, Fig. |3 we see that the allowed 
region becomes smaller and shifts towards smaller values 
of n s and fl m as A increases, in which case there is again 
a similarity with the ACDM model results. 

We have also studied the dependence of the first three 
peaks and first trough locations on parameter A. We find 
that this dependence is extremely small and can safely 
be neglected. We notice that the values of A and A 
we have previously considered all correspond to the per- 
manent acceleration regime; however, as should be ex- 
pected, changing A so as to reach the transient accelera- 
tion regime does not alter the analysis since this regime 
is not significantly different from the one where acceler- 
ation is permanent until the present, and therefore peak 
positions should not be affected. 

Hence, the model's behaviour depends essentially on 
parameter A, which measures the amount of "early 
quintessence" (using the terminology of Ref. 31]) since 
f2^ ~ 3/ A 2 , where Cl 1 ^ is the average fraction of dark 
energy before last scattering 
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n{? = r,; 1 / " fi (T)dT . (20) 
Jo 

In fact, we obtain, for A = 8, 0.039 < f^ s < 0.043, for 
the cosmological parameter range we are considering. 

We conclude that, for n s w 1, the ACDM model 
becomes increasingly disfavoured compared with the 
quintessence model of Eq. (JIJ as the amount of early 
quintessence becomes higher (A < 15). For n s < 1, the 
opposite is true i.e ACDM is favoured as compared to 
quintessence if A < 15. Notice that, as A increases (A > 
18), independently of the value of n s , the model's results 
become comparable to ACDM's, as should be expected 
since fi^ decreases. Moreover, we see that quintessence is 
distinguishible from ACDM only for h < 0.73 and n s « 1. 

Finally, we would like to mention that the non- 
negligible values we obtain for Cl s J (e.g. for A = 8, we 
get 0.022 < fiy < 0.026), typical of early quintessence 
models, will lead, as shown in Ref. [^U, to suppressed 
clustering power on small length scales as suggested by 
WMAP / CMB /large scale structure combined data. 
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APPENDIX 

We have used the analytic approximations for the 
phase shifts found in Ref. [33, which we reproduce here 
for completeness. The overall phase shift is given by 

(p = (1.466 - 0.466n s ) [ ai r° 2 + 0.291^] , (A.l) 

where 



where 

,9i = 0.0783w fc ~ a238 [l + 39.5^- 763 ]" 1 , 

.92 = 0.56[l + 21.1w b 1 - 81 r 1 . (A.5) 

The relative shift of the first acoustic peak is zero, 
6ifi = 0, and the relative shifts of the second and third 
peaks are given by 

Sif 2 = c - cm - c 2 r-° 3 + 0.05(n s - 1) , (A.6) 

with 

c = -0.1 + To.213 - 0.123fif 8 ) 

x exp j- ^52 - 63.6fif s ) u b \ , 

ci = 0.015 + 0.063 exp (-3500^) , 
c 2 = 6 x 10" 6 + 0.137(w b -0.07) 2 , 
c 3 = 0.8 + 2.3fi£ + (70 - 1260,0 , (A.7) 



6ip 3 = 10 - di< 2 + 0.08(n s - 1) , (A.8) 

with 

di = 9.97+ (3.3- 3nf s ) Lu bl 

d 2 = 0.0016 - 0.0067O* + (o.l96 - 0.22fi£) u b 

+ (2.25 + 2.770f s ) x lO -5 ^ 1 . (A.9) 

The relative shift of the first trough is 

<fy>3/2 = b Q + hrl /3 exp(6 2 r*) + 0.158(n s - 1) (A.10) 
with 



01 = 0.286 + 0.626w fc 

a 2 = 0.1786- 6.308w fc + 174.9w 2 - 1168^ (A.2) 

with u>b = 0;,/i 2 , are fitting coefficients, is given by 
Eq. (EDI) and 

r* = p r (zis) I p m (zis) (A. 3) 

is the ratio of radiation to matter at decoupling. A con- 
venient fitting formula for zi s is |48| 

z ls = 1048[1 + 0.00124 Wfc -°- 738 [l + g lW %] , (A.4) 



b = -0.086- 0.079fi£ - (2.22- 18.10^)^ 

- (140 + 4030%) to 2 , , 

bi = 0.39 - 0.98fi£ - (l8.1 - 29.2fif s ) u b 

+440w 2 , (A.ll) 
b 2 = -0.57- 3.8 exp(-2365w 2 ) , (A.12) 

Notice that the deviation of the acoustic extrema loca- 
tions calculated using the above formulae from the values 
obtained by CMBfast code is < 3% for a sufficiently wide 
range of parameters. 
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